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It is shown how to systematically construct the XX quantum spin chains with nearest-neighbor 
interactions that allow perfect state transfer (PST). Sets of orthogonal polynomials (OPs) are in 
correspondence with such systems. The key observation is that for any admissible one-excitation 
energy spectrum, the weight function of the associated OPs is uniquely prescribed. This entails 
the complete characterization of these PST models with the mirror symmetry property arising as a 
corollary. A simple and efficient algorithm to obtain the corresponding Hamiltonians is presented. 
A new model connected to a special case of the symmetric q-Racah polynomials is offered. It is 
also explained how additional models with PST can be derived from a parent system by removing 
energy levels from the one-excitation spectrum of the latter. This is achieved through Christoffel 
transformations and is also completely constructive in regards to the Hamiltonians. 



I. INTRODUCTION 

The problem of perfect state transfer (PST) in quan- 
tum information processing is deservedly attracting much 
attention. (See Refs. [5], [9] for reviews.) The transport 
of quantum state from one location to another is perfect 
if it is realized with probability 1, thereby avoiding dis- 
sipation. Occurrences of perfect transmission have been 
found in some XX chains with inhomogeneous couplings 
[5], [9], [1], [10]: in these cases, the probability for the 
transfer of a single spin excitation from one end of the 
chain to the other is indeed found to be 1 for certain 
times. These models have the advantage that the perfect 
transfer can be done without the need for active control. 
We shall here describe how such systems can be system- 
atically " constructed" from given one-excitation spectra. 

Specifically, we shall focus on Hamiltonians H of the 
XX type with nearest-neighbor interactions, i.e. 
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where J; are the constants coupling the sites I — 1 and 
I and Bi are the strengths of the magnetic field at the 
sites I (I = 0, 1, ... , iV). The symbols erf, erf , erf stand 
for the Pauli matrices which act as follows on the single 
qubit states |0) and |1): 

o*|0) = |l>, <7»|0> = -i|l>, <t z |0) = -|0) 



a*|l) = |0), ^|l) = i|0), 
It is straightforward to see that 
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and therefore that the eigenstates of H split in subspaces 
labeled by the number of spins over the chain that are 



in state |1). For our purposes, it will suffice to restrict 
H to the subspace spanned by the states witch contain 
only one excitation (or spin up). A natural basis for that 
subspace is given by the vectors 



K) = (o,o, 



i. 
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where the only " 1" occupies the n-th position. In that 
basis, the restriction J of H to the one-excitation sub- 
space is given by the following (N + 1) x (N + 1) Jacobi 
matrix 
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where it is assumed that Ji > 0. Its action on the basis 
vectors |e n ) reads 

J|e„) = J„+i|e„ + i) + B n \e n ) + J„|e„_i). (2) 

Note also that the conditions 



•/o — Jn+i — 



(3) 



are assumed. 

It has been shown [9] that the eigenvalues of J must 
satisfy a simple PST condition for the perfect state trans- 
fer to be possible. (We shall discuss this in Section 3.) 
Since they encode 3-term recurrence relations, the Ja- 
cobi matrices are diagonalized by orthogonal polynomials 
(OPs) . Spin chains allowing PST are hence in correspon- 
dence with families of OPs. 

One question we address here is the following. Given 
spectral data sets satisfying the PST condition, can the 
corresponding spin chains with the PST property be 
found? In other words, is there a procedure to explic- 
itly obtain the parameters J n and B n that determine the 



Hamiltonians? As it turns out this question has a rather 
elegant affirmative answer. Indeed, we shall see that the 
condition for a spin chain to possess the PST property is 
very simply expressed through the weight function of the 
orthogonal polynomials associated to the chain. Hence, 
given a set of eigenvalues satisfying the PST condition, 
the weight function is then prescribed. This uniquely 
fixes (up to normalization) the corresponding orthogonal 
polynomials and therefore the coefficients J„ and B n as 
their recurrence coefficients. A simple algorithm for con- 
structing spin chains with the PST property along these 
lines will be presented. 

Another constructive issue that we consider is this: can 
one obtain different spin chains with the PST property 
by performing some appropriate surgery on the spectrum 
of a model already known to generate PST? Here again 
the answer is positive and also provides tools to explicitly 
construct many spin chains with PST from a given one. 

These two natural questions have been touched upon in 
[10] and [15]. However, the link with the weight function 
of the associated OPs had so far not been stressed and 
therein lie the constructive power that our considerations 
add. 

The outline of the paper is as follows. 

In Section 2 we briefly describe standard results con- 
cerning Jacobi matrices and orthogonal polynomials. 

In Section 3, we revisit the necessary and sufficient con- 
ditions for XX spin chains to realize perfect state trans- 
fer. We derive the spectral condition already referred 
to and obtain the expression for the weight function of 
the associated polynomials on which much of the model 
characterization hinges. 

In Section 4, we offer a simple algorithm for construct- 
ing the Hamiltonians of the chains with PST from the 
spectra and the weight functions. 

In Section 5, we give two simple examples that illus- 
trate how the proposed method applies. The first is the 
well known case corresponding to a linear spectrum and 
the Krawtchouk polynomials. The second stems from 
a hyperbolic energy spectrum and leads to a novel spin 
chain Hamiltonian with PST that can be analytically de- 
scribed. 

In Section 6, we present a "surgical" procedure consist- 
ing in the removal of energy levels from the one-excitation 
spectrum of a spin chain with PST. It is shown to lead to 
other spin chains with the same property. This spectral 
surgery amounts to performing Christoffcl transforms of 
the orthogonal polynomials associated to the parent or 
original system. Since the transformed polynomials are 
explicitly known (from the general OP theory [6]), such 
is also the case for the corresponding Jacobi matrix and 
spin chain Hamiltonian. 
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II. FINITE JACOBI MATRICES AND 
ORTHOGONAL POLYNOMIALS 

Because the matrix J is Hcrmitian, there exists an 
eigenbasis \s) such that 

J\a)=x„\a), s = 0,l,...,N. (4) 

The eigenvalues x s are all real and nondegenerate x s ^ x t 
if s ^ t. 

Consider the expansion of the eigenbasis \s) in terms 
of the basis |e„) 

N 

\s) =J2 W snK). (5) 

n=0 

From (4) it is seen that the expansion coefficients W an 
satisfy 

Jn+\W s . n+ i + B n W sn + J„W s ,„-i = x s W sn . (6) 
Hence we can present them in the form 

W sn = W saX n(Xs), (7) 

where Xn(x) are polynomials satisfying the recurrence 
relation 

JnX.n — 1 

(x) (x) (8) 

and the initial conditions 

X o = l, X-i=0. (9) 
Both bases |e„) and |s) are orthonormal 

( e m| e m) — 3nm> { s \ s ) — &ss' • 

Whence the matrix W sn is orthogonal 

N 

(10) 

71=0 

and similarly 

N 

^W sn W sm = S nm . (11) 

s=0 

From (11) it follows that the polynomials Xn{x) are or- 
thonormal on the finite set of spectral points x s 

N 

^2 WsXn{x 8 )Xm{x s ) = S nm , (12) 
s=0 

where 

w s = Wl (13) 
play the role of discrete weights for the polynomials 

Xn{x). 
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Note that 



5>. = £w2, = i 

s=0 s=0 



(14) 



which amounts to the standard normalization condition 
for the weights. 

We thus have the expansions 



N 



\s) = ^ y/w^Xn(x a )\e n ) 



and similarly 



n=0 



N 



\ e n) = ^2 V™sXn(Xs)\s). 
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(16) 



s=0 



In what follows it will be convenient to work with monic 
orthogonal polynomials 



P n (x) = Jl J 2 • • • Jn Xn(x) = X n + O^"- 1 ). 

They satisfy the recurrence relation 

Pn+i(x) + B n P n {x) + U n P n -i(x) = xP n (x), 



(17) 



(18) 



where U n = J% > 0. 

The polynomials P n (x) satisfy the orthogonality rela- 
tion 



N 



/ Pn{Xs)P'm{Xs) w s — h n S n 



s=0 



(19) 



where 



U!U 2 ...U n 



Starting from the initial conditions Po = 1, P_i = 0, it is 
possible to reconstruct all monic polynomials P n {x),n = 
1, 2, . . . , N uniquely. The polynomial Pn+i(x) will be 



P N+ i(x) = (x- x Q )(x -xi)...(x- x N ). 



(20) 



From the standard theory of orthogonal polynomials [6] 
it follows that the discrete weights can be expressed as 



w. 



IN 



Pn(x s )P n+1 (x s )' 



s = 0,l,...,N. 



(21) 



In what follows we shall take the eigenvalues x s to be in 
increasing order i.e. we shall assume that 



x < x\ < x 2 < ■ ■ ■ < x N . 



(22) 



Such an ordering is always possible because, by assump- 
tion, all eigenvalues x s of the Hermitian matrix J are real 
and simple. 

One can then present the expression for P' N+1 {x s ) in 
the form 

P'n+i{x s ) = {x s - x )(x s - xi) ■ ■ ■ x (23) 
(x s - x s _i)(a; s - a; s+ i) . . . (x s - x N ) = 
(-l) N+s \P' N+1 (x s )\. 



III. NECESSARY AND SUFFICIENT 
CONDITIONS FOR PERFECT QUANTUM 
STATE TRANSFER 

A perfect quantum state transfer is realized if for some 
fixed time T > one has 



eo) 



lew), 



(24) 



where <f> is a real number. In other words, the initial 
state |eo) evolves over time T into the state |ejv) (up to 
an inessential phase factor e^). 

Expanding the states |e ) and |ej\r) in terms of the 
eigenstates \s) with the help of formula (16) we find that 
condition (24) is equivalent to 



XN(x,) = e- i +e iTx ', a = 0,l,. 



(25) 



Now, on the one hand, since the polynomial xn(x) is 
real, it follows from (25) that 



Xn{x s ) = ±1. 



(26) 



On the other hand, the orthogonal polynomials Xn{x) 
possess the interlacing property [6]. In particular, any 
zero of the polynomial xn(x) should lie between two 
neighboring zeroes x s of the polynomial xn+i(x). This 
is possible only if the values 1 and —1 alternate in (26), 
namely if 



Xn(x s ) = (-1) 



N+s 



0,1,..., N 



(27) 



where we took into account the ordering (22), formula 
(23) and the positivity of the weights w s in formula (21). 

From (25) it follows that (27) is equivalent to the con- 
dition 



f> e iTx s _ e iTrs e i4,+ir(N+2G s ) 



(28) 



where G s arc arbitrary integers. 

Whence we have the following condition for the spac- 
ings between the successive levels 



(29) 



where M a may be arbitrary positive odd numbers. 

We thus derived the necessary conditions (29) and (27) 
for a spin chain to realize perfect state transfer. 

It is easily seen that these conditions are also sufficient. 

From (27), (21) and (23) it is observed that (27) is 
equivalent to the condition 



w s 



K N 



\P' N+ l(Xs)\ 



> 



(30) 



where the normalization constant kn (chosen from the 
condition wq + w\ + ■ ■ ■ + = 1) is not essential for 
our considerations. This formula can also be found in [2] 
where it occurs in connection with the inverse spectral 
problem for persymmetric matrices. (See below.) 
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This means that only condition (29) is crucial: once 
the eigenvalues x s satisfy (29), the weights w s > are 
then uniquely determined (up to a common factor) via 
(30) and (23). In turn, the spectral points x s together 
with the weights w s are known to determine the Hcrmi- 
tian Jacobi matrix J uniquely [6]. Hence, from the sole 
knowledge of the spectral points x s (satisfying condition 
(29)) we can uniquely reconstruct the spin chain Hamil- 
tonians with the desired properties. 

IV. RECONSTRUCTION OF THE MATRIX J 
FROM SPECTRAL DATA 

Consider now the matrix J* which is obtained from J 
by a reflection with respect to the main counter-diagonal, 
i.e. 

J* = RJR, (31) 

where the matrix R (reflection matrix) is 

/ ... 1 \ 
D ... 1 

XI = 

V i o ... o o / 

The matrix J* has the same structure as the matrix J 
(i.e. J* is a Hermitian 3-diagonal matrix) 

/b* j; o \ 

Jt Bt J 2 * 0^ 
j* = J-2 B\ J3 

V . .. J* N B* N / 
with the coefficients 

B* n = B N _ n , J* = Jjv+l-n- (32) 

Consider the corresponding monic orthogonal polynomi- 
als P* (x) satisfying the recurrence relation 

P: +1 (x) + B* n P* n (x) + U* n P n ^{x) = xP* (x). (33) 

In view of (31), the matrix J* has the same eigenvalues 
x s , s = 0, 1, . . . , N as the matrix J. 

Hence the polynomials Pn(%) are orthogonal on the 
same finite set of spectral points 

N 

X] P^ s )P^{x s )w* s = h*J nm (34) 

but with another set of weights w* . There is a remarkable 
relation between the weights w s and u>* [3], [4], [16]: 

w ° w ° - Jpu^W (35) 

Now, if J is mirror-symmetric, that is J* = J, we must 
have w* — w s and in view of (35) both sets of weights 



are equal and given by (30). A Jacobi matrix J with 
mirror symmetry hence defines a spin chain with PST. 
Conversely, if the matrix J leads to a system with PST, 
the weights w s of the associated OPs must be given by 
formula (30). From (35), it then follows that w* = w s . 
This in turn implies that J* = J, in other words, that 
J is mirror-symmetric. Note that in the mathematical 
literature, the matrices with the property J* — J are 
called pcrsymmctric [2] , [7] . 

In summary, (27) is equivalent to the condition that 
the matrix J is mirror-symmetric or persymmetric. 

That (29) and the mirror-symmetry of the matrix J 
are necessary and sufficient for perfect state transfer was 
shown in [9]. 

As a result, the exploration of systems with perfect 
state transfer proceeded customarily through the search 
for Jacobi matrices with mirror symmetry. We have just 
seen that this property follows automatically from the 
prescription (30). This observation considerably simpli- 
fies the problem of finding spin chains with PST as their 
determination from the weight formula (30) avoids the 
necessity to check the mirror symmetry of J. 

At this point, to complete the practical recipe for the 
construction of spin chains with PST, we only need to 
provide an efficient algorithm to determine the matrix J 
from the spectrum. 

One such possible algorithm (to reconstruct persym- 
metric matrices from spectral data) was given in [2] . We 
here propose another one which seems more efficient es- 
pecially in view of the "spectral surgery" procedure to 
be presented in the next section. 

First, we notice that the polynomial xn{x) can be re- 
constructed from (27) by the ordinary Lagrange interpo- 
lation procedure [8]. 

Indeed, xn{%) has degree N and takes N + 1 pre- 
scribed values ±1 at N + 1 prescribed distinct points 

Xq,Xi, . . . ,X N . 

It follows that 

N 

XN ( X ) = J2(-1) N+S £s, (36) 
where C s are the standard Lagrange polynomials 

t=0 Xs Xl 

(as usual, the symbol ' means that i ^ s in the product 
(37)). The monic polynomial P/v(x) is then obtained 
through the division of xn(x) by the coefficient of its 
leading monomial. 

We thus know explicitly two monic polynomials: 
Pn(x) and Pn+i(x) — (x — xq)(x — x\) . . . (x — Xn)- 
Starting from these polynomials it is possible to recon- 
struct step-by-step all the orthogonal monic polynomials 
P n (x), n = N - 1, N - 2, . . . , 1 by the well known Eu- 
clidean algorithm. 



■5 



Let us divide the polynomial Pn+i{x) by the polyno- 
mial Pn(x): 

P n+ i(x) = q N {x)PN{x) + Rn-i(x), (38) 

where q n (x) = x — (3n and Rn-i(x) is the residue of 
this division. By construction, we have on the one hand 
deg(i?jv_i(x)) = N — 1 and hence 

Rn-i(x) = inQn-i(x), (39) 

where Qn-\{x) = x N ~ 1 +0(x N ~ 2 ) is a monic polynomial 
of degree N — 1 . 

On the other hand, we have the recurrence relation 
(18) from which we conclude that /3/v = 6jy, In = un 
and Qn-i{x) = Pn-i(x). We thus get the next monic 
orthogonal polynomial Pn-i(x) as well as the recurrence 
coefficients b Nl u N . 

The same steps can then be repeated with the polyno- 
mials Pjv(x), Pn-i(x) yielding the recurrence coefficients 
&jv-i)^Ar-i an d the polynomial Ppj_2(x). Iteration will 
provide all recurrence coefficients b n , n = 0, 1, . . . , N and 
u n , n — 1, 2, . . . , N together with the corresponding or- 
thogonal polynomials P n (x), n = 1, 2, . . . , N. 

From a physical point of view, the case B n — is of 
special interest because it corresponds to zero external 
magnetic field. Jacobi matrices J with zero diagonal en- 
tries B n — correspond to symmetric orthogonal poly- 
nomials satisfying the property [6] 

Pn(-x) = (-l) n P n (x). (40) 

Conversely, (40) is equivalent to the condition B n = in 
the Jacobi matrix J [6]. 

What is more important is that (40) is tantamount to 
the following spectral properties [6]: 

(i) the eigenvalues are anti-symmetric, i.e. 

x n = -x N - n , n = 0,l,...,N (41) 

and 

(ii) the weights are symmetric: 

w s = w N - s > 0. (42) 

We can now apply these observations to the PST prob- 
lem. Assume that the eigenvalues x s satisfy the prop- 
erties (29) and (41). Then, the weights w s constructed 
from formula (30) obviously satisfy (42). 

For perfect state transfer to be achieved in spin chains 
with zero magnetic field it is thus necessary and sufficient 
that the one-excitation energies satisfy (29) and (41). A 
similar result was obtained in [9] using a method asso- 
ciated to the inverse eigenvalue problem. Note that it 
is assumed in [9] that the matrix J is mirror-symmetric. 
We already observed that this assumption is superfluous 
as the mirror symmetry of J follow from (30). 



V. TWO EXPLICIT EXAMPLES 

As illustrations of how the constructive procedure we 
have described can be applied, we present in this sec- 
tion two examples. The first example is well known, the 
second one seems to be new. 

Before proceeding, note that given a spectral data set 
x s , s = 0, 1, . . . , N it is straightforward to obtain another 
set satisfying (29) by an affine transformation 

x s = ax s + 0, (43) 

where a, ft are arbitrary real parameters. The corre- 
sponding Jacobi matrix will then have for its entries 

B n = aB n + p, J n = aJ n . (44) 

We can use this freedom to choose the most convenient 
form of the spectral data. 

In particular, it is always possible to choose the pa- 
rameters a, f3 so that the eigenvalues x s are integers with 
alternating parity (all X2 S even and all X2 S +i odd). We 
will use this observation in the following. 

Consider first the uniform grid 

x s = s-N/2, s = 0,l,...N. (45) 

Using (30) we easily reconstruct the weights w s to obtain 
the binomial distribution. The corresponding orthogo- 
nal polynomials P n (x; N) are the symmetric Krawtchouk 
polynomials and the entries of the matrix J are 

B „=0, J ; = "( JV+ 4 1 -' . (46) 

This example is well known and was in fact obtained in 
[1] as a first example of inhomogeneous spin chain with 
the perfect state transfer property. 

The next example is less trivial and seems to have es- 
caped notice. Take a "hyperbolic" analogue of the uni- 
form spectrum, i.e. 

x s = A(q- S+N / 2 - q- N ' 2+s ), s = 0,l,...,N, (47) 

where < q < 1 and A are real parameters. These pa- 
rameters should be chosen so as to enforce the condition 
that all differences x s+ i — x s are positive odd integers. 
To that end, it is sufficient to demand that all x s be in- 
tegers with alternating parity, say xq is even, x\ odd, x-i 
even and so on. 

To do this, let us notice that the spectral points (47) 
satisfy the recurrence relation 

x s+1 +x s _i = {q + q-^Xs, s = l,2,...,JV-l. (48) 

In order to ensure that all integers with alternating 

parity we need to require 

q + q-*=K, (49) 

where K = 4, 6, . . . is an arbitrary positive even integer 
in the case of even N and K = 6, 10, 14, 18, . . . in the case 
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of odd N (we avoid the case K — 2 because it leads to 
the degenerate case q = 1 corresponding to the uniform 
grid (45)). Condition (49) means that q is a special case 
of quadratic irrationality. The difference in the ranges of 
K for N even and N odd is explained by the observa- 
tion that for N even the minimal distance between the 
eigenvalues x s is A{q~ x — q) while for N odd the minimal 
distance is 2A(q~ x / 2 — q 1 / 2 ). The admissible values for 
K then easily follow from (48). 

With the help of (30) the weights are straightforwardly 
found (up to a normalization factor): 



w s = (-l) s q s "(l + q 2s -») 



(q;q) s (-q;q) s 



(50) 



where 



(x; q) n = (1 - x)(l -xq)...(l- xq r ' 



stands for the q-shifted factorial [11]. 

The weights (50) correspond to a special case of the q- 
Racah polynomials that are orthogonal on the grid (47) 
[11]. Hence the entries B n ,J n of the Jacobi matrix J 
can be recovered from the known recurrence coefficients 
of the q-Racah polynomials [11]. 

The spectrum x s is anti-symmetric xn-s = —x Sl hence 
the diagonal terms are absent B n = 0. For the non- 
diagonal terms we have from [11] the expression 



U n = Jl = A 2 (l-*? 2 ")^"-^-!) 



(l + q 



2n-N-2 



){1+q 2n- N y 



It is easily seen that the coefficients (51) are positive and 
satisfy the mirror symmetry condition U n — Un+i-u- 



VI. SPECTRAL SURGERY AND GENERATING 
NEW PST CHAIN 

Given a one-excitation spectrum x s , s = 0,1,..., N 
satisfying condition (29), we have seen how to construct 
the spin chain Hamiltonian with PST. The Jacobi matrix 
which has the couplings and magnetic field strengths as 
its entries is recovered from the prescribed weights w s 
given by (30), or, in an equivalent form 



{-l) N+s 



11=0 ( X s 



Xi) 



(52) 



where it is assumed that x < x\ < ■ ■ ■ < x^. By con- 
struction, all weights are positive w s > 0, s = 0, 1, . . . , N. 
The weights w s in (52) are defined up to an arbitrary pos- 
itive common factor which has no effect on the entries 
Ji, Bi of the Jacobi matrix J. 

Consider the modified set of spectral data 
xi, X2, ■ ■ ■ , xn obtained by removing the first entry 
xq. The corresponding weights 



w s = 



(-1) 



N+s 



( x s x i) 



8 = 1,2,. ..,N 



can be obtained from the initial weights in the following 
simple manner: 



w s = const(x s — xo)w s . 



(53) 



This procedure removes the eigenvalue x$ and preserves 
the positivity of the resulting weights w s . 

Similarly, one can remove any fixed spectral point xf 

w s = const(x s — Xj)w s , s = 0, 1, . . . , j — l,j + 1, . . . , N. 

(54) 

In this case, however, the new weights will be positive 
only if either j = or j = N. In all other cases the 
weights w s cannot be positive for all s. 

Nevertheless, removing a pair of neighboring points 
Xj,Xj + i through 



w s = const(x s — Xj)(x s — Xj + i)w s . 



(55) 



maintains positivity of the weights w s for all s, preserves 
property (29) and thus provides a new admissible Jacobi 
matrix J which also generates perfect state transfer. 

This procedure of removing pairs of neighboring levels 
can obviously be iterated. 

The orthogonal polynomials P n (x) corresponding to 
the weights (54) are obtained from the polynomials P n (x) 
by the Christoffel transform [14] 



Pn{x) 



Pn+l(x) - A n P n (x) 



(56) 



(51) 

where 



A n — 



n+l 



Pn (Xj ) 

The entries of the matrix J are related to the entries of 
the matrix J by the well known formulas [18] 

A n 



U n = Un 



A n ~ 



Bn 



B n +i + A n+ i — A n . (57) 



Formulas (56) and (57) can be applied iteratively in order 
to obtain new matrices J with perfect state transfer from 
a matrix J with that property. 

Assume now that the initial matrix J describes a dy- 
namics with zero magnetic field, i.e. B n = 0. The spec- 
tral points x s satisfy therefore the symmetry condition 
(41). In this case one may remove a symmetric pair of 
boundary eigenvalues corresponding to the first and last 
levels 



w s = const (x 8 — Xq)w s 

and find for the associated polynomials P n (x) [14] 

-^1+2(^0) 
P n (xo) 



Pn{x) = 



Pn+2(x) - K n P n (x) 



K n = 



(58) 



(59) 



therefore obtaining a new Jacobi matrix J with zero mag- 
netic field B n = and with 



U n = U n 



K, 



n-l 



(60) 
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If N is odd, there is another possibility to remove two 
neighboring levels, in this case from the middle of the 
spectrum. Indeed, take 

w s = const(x s — Xj)(x s — Xj + i)w s = const{x 2 s — Xj)w s , 

(61) 

where j = (N — l)/2. Formulas (59) and (60) remain 
valid if one replace x with Xj . 

The idea of obtaining PST chains with B n = from 
given ones was discussed in [15]. Our approach is much 
more explicit and exploits well known formulas from the 
theory of orthogonal polynomials. 

As a simple (but already nontrivial) example, consider 
the removal of levels from the middle of the uniform spec- 
trum (45) with odd N. After a finite number L of such 
transformations, we obtain spectral data of the form 

x s = -N/2, -N/2 + 1, L — 3/2, —L - 1/2, 

L + 1/2, L + 3/2,..., N/2 - 1, N/2. (62) 

For L = this spectral data set coincides with (45). For 
L = 1,2, ■■■ < (N — l)/2 we have two uniform grids 



separated by a gap of length 2L + 1. This spectrum 
corresponds to a spin chain with B n = that has the 
perfect transfer of state property and that was introduced 
in [12] and studied in [13]. We have shown in [17] that 
this model is associated to the dual -1 Hahn polynomials. 

Let us conclude by stressing here that ALL spin chains 
with the PST property (not only those with B n = as 
in [15]) can be obtained through such (iterated) spectral 
surgery from the PST spin chain corresponding to the 
uniform grid (45). 
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